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ABSTRACT 

Most carbon-enhanced metal-poor (CEMP) stars are thought to result from past mass transfer of He-burning material from an asymp- 
totic giant branch (AGB) star to a low-mass companion star, which we now observe as a CEMP star Because AGB stars of intermediate 
mass efficiently cycle carbon into nitrogen in their envelopes, the same evolution scenario predicts the existence of a population of 
nitrogen-enhanced metal-poor (NEMP) stars, with [N/Fe] > 1 and [N/C] > 0.5. Such NEMP stars are rare, although their occurrence 
depends on metallicity: they appear to be more common at [Fe/H] < -2.8 by about a factor of 10 compared to less metal-poor stars. 
We analyse the observed sample of metal-poor stars with measurements of both carbon and nitrogen to derive firm constraints on the 
occurrence of NEMP stars as a function of metallicity. We compare these constraints to binary population synthesis calculations in 
which we vary the initial distributions of mass, mass ratio and binary orbital periods. We show that the observed paucity of NEMP 
stars at [Fe/H] > -2.8 does not allow for large modifications in the initial mass function, as have been suggested in the literature to 
account for the high frequency of CEMP stars. The situation at lower metallicity is less clear, and we do not currently have stellar 
models to perform this comparison for [Fe/H] < -2.8. However, unless intermediate-mass AGB stars behave very differently at such 
low metallicity, the observed NEMP frequency at [Fe/H] < -2.8 appears incompatible with the top-heavy forms of the initial mass 
function suggested in the literature. 

Key words, stars: carbon - binaries: close - stars: chemically peculiar - Galaxy: halo - Galaxy: stellar content - nucleosynthesis, 
abundances 



1. Introduction 

Carbon-enhanced metal-poor (CEMP) stars make up a large 
proportion of the most metal-poor stars in the Galactic halo. 
Estimates of the fraction of very metal-poor stars (with [Fe/H] < 
-2) having carbon enhancements [C/Fe] > 1.0 range between 
10% and 25 % ( |Frebel et al.|2006l [Lucatello et al.|2006||Carollo| 
|et al. 20121. Based on their abundance patterns different sub- 



groups are identified among the CEMP stars (Beers & Christlieb 
|2005 [Masseron et al.|2010j . By far the most numerous group 
are also enriched in s-process elements, which make up about 
80% of CEMP stars ( |Aoki et al.||2007l ). We refer to this group 



as CEMP-s stars, noting that a certain fraction of these are also 
enriched in r-process elements (these have been designated as 
CEMP-rs stars in the literature). Radial velocity monitoring sug- 
gests that probably all CEMP-s stars are binaries ( [Lucatello et al.| 



2005b I. A widely accepted formation scenario for this group in- 



volves pollution by mass transfer in a binary system from a more 
massive asymptotic giant branch (AGB) companion, which has 
since become a white dwarf. Of the CEMP stars without strong 
s-process enhancements, one star is found to be strongly en- 
riched in r-process elements, the remainder having no or only 
weak overabundances of heavy elements. The latter group are 
known as CEMP-no stars. The origin of these stars is less clear; 
there is no strong evidence for a large binary fraction among the 



CEMP-no stars. They appear to be more common at the lowest 
metallicities. Although they may have followed a similar evolu- 
tion path as the CEMP-s stars, it is also possible that they have 
an entirely different origin (e.g. Meynet et al.|2010] l. 



Within the mass transfer scenario, the large proportion of 
CEMP-s stars requires the existence of a sufficient number 
of binary systems with primary components that have under- 
gone AGB nucleosynthesis. Several authors ( Abia et aL][2001 



[Lucatello et al.,2005a^ Komiya et al.[2007) have therefore argued 
that a different initial mass function (IMF), weighted towards 
intermediate-mass stars, is needed at low metallicity. Lucatello 
et al. (2005a ) find that the standard IMF for the Galactic disk 
( [Miller & Scalo| 1979| l contains too few AGB progenitors in 
the mass range 1.2-6 Mq, assumed to be responsible for car- 
bon and .s-element production, to explain the observed fraction 
of CEMP-s stars. They propose an IMF with a mild shift to- 
wards intermediate-mass stars is required. Komiya et al. (2007i 
present a semi-analytical model for binary population synthesis, 
based on AGB models in which low-mass stars (0.8-1.5 Mg) at 
[Fe/H] < -2.5 produce carbon and s-elements by dual shell 
flashes (see Sect. |5]for details). While this increases the param- 
eter space for making CEMP-s stars considerably, it is barely 
sufficient to produce a CEMP-s fraction of more than 10% with 
a Galactic -disk IMF. They further propose that the CEMP-no 
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Stars are produced in binary systems with AGB primaries of 
3.5-6 Mq. From an assumed CEMP-no/CEMP-s ratio of at least 
1/3, consistent with observations available to them, they find 
consistency with both the CEMP-s fraction (10-25%) and the 
CEMP-no/CEMP-s ratio for a top-heavy IMF, with a median 
mass of 10 Mq. Such proposed changes to the IMF would have 
profound consequences, since the IMF affects such diverse quan- 
tities as the mass-to-light ratio of galaxies, the overall production 
of heavy elements and galactic chemical evolution. 

However, the model calculations on which these conclusions 
are based still contain many uncertainties regarding the evolu- 
tion and nucleosynthesis of low-metallicity AGB stars, the effi- 
ciency of mass transfer, and the evolution of the surface abun- 
dances of the CEMP stars themselves. A different approach to 



the problem of ubiquity of CEMP stars was taken by Izzard 



et al. (2009 1. They performed detailed binary population synthe- 



sis ca culations to explore whether the observed CEMP-s frac- 
tion at [Fe/H] ^ -2.3 can be explained within the boundaries 
of the model uncertainties, while assuming a standard solar- 
neighbourhood IMF. Their results confirmed that with current 
detailed AGB models, in which dredge-up of carbon occurs only 
for masses larger than 1 .2 Mq, the large CEMP fraction cannot 
be explained. However, they also found that efficient dredge- 
up in lower-mass AGB stars (0.8-1 .2 Mq) can substantially in- 
crease the expected CEMP fraction to within the observed range. 
Therefore the question of whether CEMP stars provide evidence 
for a modified IMF at low metallicity remains open. 

Apart from carbon, substantial enhancements of nitrogen 
with respect to iron are common among CEMP stars, typically 
with [N/C] < 0. Detailed AGB nucleosynthesis models of low 
initial mass (< 3-5 Mq, depending on metallicity) make carbon 
through the dredge-up of helium-burning products, but produce 
very little nitrogen. On the other hand, AGB stars of higher mass 
are predicted to convert the dredged-up carbon into nitrogen by 
CN cycling at the bottom of the convective envelope (hot bottom 
burning, HBB). The surface abundances of these more massive 
AGB stars approach the CN equilibrium ratio of [N/C] 2. 
Detailed evolution models of AGB stars ( |Karakas & Lattanzio] 
2007 1 indicate that, while at solar metallicity HBB sets in at 
^ 5 Mq, significantly lower masses are needed at low metal- 
licity. At [Fe/H] - -2.3 the lower mass limit is found to be 



ss 3Mq ( Karakas 2010 1. By the same mass-transfer scenario 
that forms CEMP-s stars, one may thus expect a population of 
so-called nitrogen-enhanced metal-poor (NEMP) stars, with en- 
hanced nitrogen abundances and [N/C] > 0.5. Note that we thus 
expect NEMP stars to result from a similar AGB mass range 
that Komiya et al. ( 2QQ7\ suggested to produce CEMP-no stars 
(based on models in which no HBB occurs). Although a few ex- 
amples of NEMP stars are known, mostly at [Fe/H] < -2.9, they 
appear to be very rare ( Johnson et al.|200T i. 

In this paper we show that the number of NEMP stars sets 
an additional constraint on possible changes to the IMF at low 
metallicity. In Sect. 2 we discuss the observational evidence re- 
garding NEMP stars. Our population synthesis models are de- 
scribed in Sect. 3 and the results are presented in Sect. 4. We 
discuss our results in Sect. 5 and in Sect. 6 we give our conclu- 



2. Nitrogen-enhanced metal-poor stars 



I Johnson et al.| ( |2007| l used the same theoretical considerations 
as we gave above to argue for the expected existence of a sub- 
stantial population of NEMP stars, which they defined as stars 
with [N/Fe] > 0.5 and [N/C] > 0.5. They set out to search 
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Fig. 1. Distribution of [N/Fe] versus [C/Fe] for metal-poor stars 
with observed carbon and nitrogen abundances. Star-shaped 
(red) symbols are stars with [Fe/H] < -2.8, triangles (dark 
blue) have -2.8 < [Fe/H] < -1.8 and circles (yellow) have 
[Fe/H] > -1.8. Giants (stars with logg < 4) are shown as filled 
symbols while higher-gravity stars have open symbols. Large 
symbols represent data from Masseron et al.| (2010) and the 
SAGA database; data from |Lucatello et al.| ( |2006 1 is shown with 
smaller symbols. Multiple data for the same star are connected 
with solid line segments. The dashed and solid lines indicate the 
criteria adopted for CEMP stars and NEMP stars, and the dash- 
dotted line is the alternative criterion [(C + N)/Fe] > 0.9 (see 
text). 



for NEMP stars by determining nitrogen and carbon abundances 
for a sample of 21 mildly carbon-enhanced stars. They found 
no NEMP stars but confirmed the known trend that [N/C] val- 
ues of CEMP stars tend to be intermediate between what is ex- 
pected from AGB stars with masses above and below the hot 
bottom burning limit. While their study showed that NEMP stars 
are rare, and their sample has the advantage of being homoge- 
neous, it is too small to derive firm constraints on the frequency 
of NEMP stars. We therefore performed a literature search to 
obtain a more complete overview of the situation. 

Characterising a star as NEMP requires measurements of 
both the C and N abundances, which are only available for a 
limited number of halo stars. In Fig. [T] we show the distribution 
of [N/Fe] versus [C/Fe] of 259 metal-poor stars with measure- 
ments of both carbon and nitrogen. Data for 176 stars have been 
col lected in the SA GA database ( |Suda et al.||2008l [2011] ) and 
by Masseron et al. ( 2010| ). These datasets largely overlap but 
Masseron et al.| list~a few stars not included in SAGA. Many stars 
have multiple measurements of both C and N that are broadly 
consistent within the observational errors, in which case we gen- 
erally take the data from the most recent publication as listed by 
Massero n et al.|p010[ ). For a few stars with clearly discrepant 
measurements, we show two data points connected by a line 
segment. Added to this selection are 87 stars, mostly CEMP, 
with C and N abundances determined by |Lucatello et al.| ( |2006| l. 
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Table 1. Nitrogen-enhanced metal-poor stars, and other possible remnants of intermediate-mass binaries (last two entries). 



star 


[Fe/H] 


logs 


[C/Fe] 


[N/Fe] 


L/C 


[Na/Fe] 


[Mg/Fe] 


[Ba/Fe] 


reference 


CS22949-037 


-3.86 


1.7 


+ 1.07 


+2.11 




+ 1.45 


+ 1.50 


-0.59 


COS 




-3.97 


1.5 


+ 1.17 


+2.57 


3 


+2.09 


+ 1.58 


-0.58 


D02 


CS22960-053 


-3.14 


2.1 


+2.05 


+3.06 






+0.65 


+0.86 


A07 




-3.08 


2.4 


+ 1.15 


+ 1.15 










J07 


CS29528-041 


-3.30 


4.0 


+ 1.61 


+3.09 




+ 1.20 


+0.45 


+0.97 


S06 


CS303 14-067 


-2.85 


0.7 


+0.50 


+ 1.20 


5 


-0.08 


+0.42 


-0.57 


A02 




-2.76 


1.0 


+0.25 


+0.50 










J07 


CS30322-023 


-3.39 


-0.3 


+0.80 


+2.91 


4 


+ 1.29 


+0.80 


+0.52 


M06, MIO 




-3.25 


1.0 


+0.56 


+2.47 




+ 1.04 


+0.54 


+0.59 


A07 


HD25329 


-1.8 


4.7 


+0.11 


+1.00 


>40 


+0.24 


+0.59 


+0.11 


GOO, G03, FOO 


HD206983 


-0.99 


0.6 


+0.50 


+ 1.21 


5 




+0.26 


+0.92 


MIO 


HE0400-2030 


-1.73 


3.5 


+ 1.14 


+2.75 




+0.71 


+0.62 


+ 1.64 


A07 




-1.7 


2.5 


+0.8 


+1.0 










L06 


HE 103 1-0020 


-2.93 


2.2 


+ 1.73 


+2.43 






+0.53 


+ 1.28 


C06 


HE1337+0012" 


-3.20 


4.4 


+0.49 


+ 1.42 




-1.13 


+0.41 


-0.25 


A06 


HE1410+0213 


-2.52 


2.0 


+2.33 


+2.94 


3 




+0.33 


+0.05 


MIO 




-2.23 


3.5 


+ 1.83 


+ 1.73 






+0.21 


+0.14 


C06 


HE1413-1954 


-2.9 


4.0 


+ 1.7 


+2.5 










L06 


HE2 150-0825 


-1.3 


3.8 


+0.3 


+ 1.3 










L06 


HE2253-4217 


-2.4 


1.6 


+ 1.0 


+ 1.5 










L06 


CS29528-028 


-2.86 


4.0 


+2.77 


<+3.58 




+2.68 


+ 1.69 


+3.27 


A07 


SDSS 1707+58 


-2.52 


4.2 


+2.1 






+2.71 


+ 1.13 


+3.40 


AOS 



References. A02:|Aoki et al. (2002), A0 6:|Aoki et"ari([2006ll, A07:|Aoki et allpOOTll, A08:|Aoki et al lpOOSll, C06:|Cohen et al |(|2006l>, COS 



[Cohen et al.|(l2008|l, D02r|Depagne eTai] ( |2002|l, F00:|Fulbright](j2000| l, G00:|Gratton et aTjpOOO) , G03:|Gratton et al.|i|2003| l, J07: |johnson etaT] 
\2001) , L06: |Lucatello et al.H2006| l, M06: |Masseron et al.H2006| , M10: |Masseron et al.| ( |2010t , S06: |Sivarani et al.| ( |2006t . 

^^'^so known as G64-12. 



which were not included in the SAGA database. We note that the 
[Johnson et'aT c'IOO?) data are included in SAGA. 

Following Izzard et al. ( 2009 hereinafter Paper I), we desig- 
nate as CEMP stars those with [C/Fe] > 1.0 and as NEMP stars 
those with 



[N/Fe] > 1.0 and [N/C] > 0.5. 



(1) 



Note that these definitions partly overlap; stars that satisfy both 
criteria (i.e. stars with [C/Fe] > 1.0 and [N/C] > 0.5) are desig- 
nated as CNEMP stars. The adopted CEMP and NEMP criteria 
are shown as dashed and solid lines in Fig.[T] 

Our NEMP criterion ([T]) is somewhat more restrictive on the 
nitrogen abundance than the definition of [Johnson et al.[ (2007 1. 
Application of this criterion to our selection of datavields 24 
NEMP star candidates (see also Table 2 of Paper v). Ten of 
these stars, however, have [C/Fe] ^ 0.0 and [N/Fe] values be- 
tween 1.0 and 1.3. All these stars have [Fe/H] < -2.8 and are 
on the upper part of the red giant branch (RGB): they belong to 
the sample of 'mixed stars' of Spite et al. (2005 2006 1. Their 



nitrogen enhancements are very likely the result of CN cycling 
and 'extra mixing' on the RGB itself, and not the result of mass 
accretion of nitrogen-iich material. Indeed, their combined C + 
N abundances are only mildly super-solar, [(C + N)/Fe] $ 0.7. 
Our aim of eliminating such stars from our selection by the cri- 
terion [N/Fe] > 1.0 was obviously not completely successful. 

An alternative, and perhaps better, criterion to distinguish 
stars that have had significant accretion of He-burning products 
(either carbon or nitrogen) can be based on the C + N abundance, 
because CN-cycling during HBB conserves the total number of 
C + N nuclei. If carbon is enhanced by a factor greater than 10 



' We found 17 stars in Paper I. The difference is due partly to the 
addition of new data and partly to our use in Paper I of the average 
observed abundance, while here we use data from a single source. 



(corresponding to our CEMP criterion [C/Fe] > 1.0) and sub- 
sequently cycled partly into nitrogen, we expect [(C + N)/Fe] > 
0.9, assuming [Grevesse et al.[ ( |I996| l solar-scaled initial abun- 
dances. This criterion is shown as the dash-dotted line in Fig.[T 
It effectively eliminates stars that are N-rich owing to interna 
mixing processes on the RGB. In this paper, however, we pre- 
fer to use our criterion ([T]i because it is consistent with Paper I 
and our model calculations, but we will subject NEMP stars with 
[(C + N)/Fe] < 0.9 to critical examination. 



If we eliminate the 10 'mixed' RGB stars we are left with 
14 NEMP candidates, for which we list several key properties 
in Table[T] However, not all of these stars necessarily qualify as 
'genuine' NEMP stars, i.e., stars that are likely to have been pol- 
luted by nitrogen-rich material, possibly from an intermediate- 
mass AGB star Other clues can be obtained from the '^C/'^C 
ratio, which is expected to be close to its CN-equilibrium value 
of 3-4 for material that has gone through HBB, and from mea- 
surements of Na and Mg. In an AGB star undergoing HBB, 
Na is produced by the NeNa chain from the abundant ^^Ne 
brought to the surface by third dredge-up, along with Mg if the 
-^Ne(Q',n)^^Mg source is active. Therefore [Na/Fe] and [Mg/Fe] 
values greater than zero are expected in material lost by an 
intermediate-mass AGB star (e.g. see Lugaro et al. 20I2[ l, al- 
though we note that Na production by AGB stars is uncertain 
owing to the large uncertainty in the ^^Ne(p,y)^^Na reaction rate 
( [Izzard et al."|2007) . Clearly, not all NEMP candidates in Table [T] 
for which ' C/'^C, Na and/or Mg are measured show the ex- 
pected behaviour. In particular, HD25329 is an unevolved main- 
sequence star with [N/Fe] = 1.0 and [C/Fe] = 0.1, and thus 
a modest [(C + N)/Fe] = 0.5. Its '^c/'^c ratio is at least 40 
which is inconsistent with accretion of material that has under- 
gone HBB. Although the origin of its high N abundance is un- 
clear, we should probably exclude this star from the sample of 
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Table 2. Incidence of CEMP, NEMP and CNEMP stars among stars with carbon and nitrogen measurements. The columns labelled 
/max give the 90%, 99% and 99.9% upper limits to /nemp and /cnemp for the selections with \ogg < 4. 



selection 






CEMP 


NEMP 


CNEMP 


/nemp 


/max(90%/99%/99.9%) 


/cnemp 


/max(90%/99%/99.9%) 


all 






156 


14 


8 


0.086 




0.051 




all, log g < 4 






142 


12 


8 


0.082 


0.122/0.156/0.185 


0.056 


0.092/0.123/0.149 


[Fe/H] < -2.8 






26 


8 


5 


0.276 




0.192 




[Fe/H] < -2.8, log 


8 


< 4 


24 


7 


5 


0.269 


0.411/0.514/0.588 


0.208 


0.352/0.459/0.539 


-2.8 < [Fe/H] < - 


1. 


8 


101 


3 


2 


0.029 




0.020 




-2.8 < [Fe/H] < - 


1. 


8,logg<4 


92 


2 


2 


0.022 


0.057/0.088/0.116 


0.022 


0.057/0.088/0.116 


[Fe/H] > -1.8 






29 


3 


1 


0.097 




0.034 




[Fe/H] > -1.8, log 


8 


< 4 


26 


3 


1 


0.107 


0.223/0.316/0.390 


0.038 


0.142/0.229/0.304 


[Fe/H] > -2.8 






130 


6 


3 


0.045 




0.023 




[Fe/H] > -2.8, log 


8 


< 4 


118 


5 


3 


0.042 


0.076/0.106/0.131 


0.025 


0.056/0.083/0.106 



genuine NEMP stars. We comment in more detail on the stars in 
Table [T] in Appendix [A] 

Fig. [T] demonstrates the paucity of NEMP stars, in particu- 
lar among the more metal-rich part of the population. We focus 
in particular on the metallicity range -2.8 < [Fe/H] < -1.8 
because these stars can be best compared to our models at 
[Fe/H] = -2.3. Within this metallicity range (triangles) there 
are only three candidate NEMP stars, one of which (HD25329) 
we can most likely exclude. The other two are also doubtful 
cases: HE1410-I-0213 has an alternative measurement that gives 
[N/C] < 0, and HE2253-4217 is exactly on the borderline of our 
criterion with [N/C] = 0.5 and [C/Fe] = 1.0. Thus, there may be 
no NEMP stars in the [Fe/H] range that we can directly compare 
with our models at [Fe/H] - -2.3. In contrast, NEMP stars ap- 
pear to be much more common among the low-metallicity pop- 
ulation ([Fe/H] < -2.8, star symbols). However, we are dealing 
with low-numbers statistics so we need to take care in drawing 
such conclusions. We therefore perform a more careful analysis 
of the statistics of NEMP stars. 



2. 1. Statistics of CEMP and NEMP stars 

It is difficult to derive an overall fraction of NEMP stars from the 
observed samples. The latest release of the SAGA database of 
metal-poor stars (dated Oct. 31, 2010; |Suda et al.|2011 



contams 

656 stars with measurements of the carbon abundance, of which 
1 10 (about 17%) classify as CEMP stars. For 148 stars the niti-o- 
gen abundance has also been determined. Out of these, 68 have 
[C/Fe] > 1, amounting to a CEMP fraction of 46% among stars 
with both C and N measurements. Apparently there is a strong 
bias towards selecting CEMP stars as targets for nitrogen abun- 
dance determinations. The Masseron et al. (2010) and LucatellOi 



et al. ( 2006 1 samples of stars with C and N measurements are 
even more strongly biased towards CEMP stars. One reason for 
this bias may be that in most studies nitrogen abundances are 
determined from the CN band, which requires the presence of a 
substantial amount of carbon. 

On the other hand, among stars for which both C and N 
abundances are measured there is no reason to suspect a strong 
bias towards certain [N/C] values. We can thus expect that the 
number ratio of NEMP to CEMP stars among these stars is a 
fairly unbiased observable. We designate as /nemp the fraction 
of NEIVIP stars among stars that ai-e either CEMP or NEMP 
(i.e. /nemp = A^NEMp/[A^CEMP + A^NEMP - A^CNEMp], which avoids 
counting CNEMP stars twice) such that /nemp < 1 by definition. 
The common use of CN as a nitrogen abundance indicator may 
bias the observed sample of NEMP stars towards C-rich stars in 



such a way that non-C-rich NEMP stars are underrepresented in 
the observed sample. We therefore also consider the fraction of 
CNEMP stars among CEMP stars, which we denote as /cnemp 
and which should be unaffected by such a bias. Another advan- 
tage of using the (C)NEMP/CEMP ratio as a constraint on our 
model predictions is that it is independent of various model un- 
certainties, most notably the binary fraction. We therefore at- 
tempt to derive /nemp and /cnemp from the 162 stars that are 
classified as either CEMP or NEMP from the combined SAGA- 
Masseron-Lucatello sample as shown in Fig. [T] From the 14 
NEMP stars and 8 CNEMP stars we have found in this sam- 
ple (Table [TJ we thus deduce a nominal /nemp = 0.086 and 
/cnemp = 0.05 1 in the entire sample. 

However, this number deserves closer scrutiny for two rea- 
sons. Firstly, Fig. [T] and Table [T] demonstrate that most NEMP 
stars have rather low [Fe/H] and the NEMP/CEMP ratio appears 
to depend on metallicity. This dependence becomes clear when 
we split up the sample into different [Fe/H] ranges, see Table [2] 
Among stars with [Fe/H] < -2.8 both /nemp and /cnemp are 
about 10 times higher than in the range -2.8 < [Fe/H] < -1.8. 

Secondly, the numbers quoted in Table |2] include stars of un- 
certain NEMP status, so we can at best derive upper limits to the 
NEMP/CEMP fraction from the observations. Since we are deal- 
ing with very small numbers, especially when we consider sub- 
samples in different metallicity ranges, we cannot simply take 
the nominal fractions in Table|2] Instead we should ask the ques- 
tion: given a true fraction / of NEMP stars among CEMP stars, 
what is the probability of finding at most NEMP stars in an ob- 
served sample of A^tot CEMP stars? Assuming a purely stochastic 
process, this probability is P{< N) = Yjn=o ^bin,Ntot,f), where 
Pb(n,Ntoi,f) is the binomial probability of drawing n NEMP 
stars out of a sample of A^tot given a true fraction /. Since P(< N) 
is a decreasing function of /, if we find P{< N) to be, say, 10% 
for a certain value of /, we can exclude NEMP/CEMP fractions 
larger than / with 90 per cent confidence. In this way we can 
compute upper limits to / for certain confidence limits, i.e. 90%, 
99% and 99.9%. These upper limits are given in Table|2]for the 
different subsamples and for both /nemp and /cnemp- 



As in |Izzard et al. (2009 1 we limit our analysis to stars with 
logg < 4, i.e. subgiants and giants. This facilitates the com- 
parison with our model calculations and does not introduce sig- 
nificant biases in the numbers we derive, as can be verified in 
Table |2] We further limit the selection to stars with -2.8 < 
[Fe/H] < -1.8, which is the relevant metallicity range for com- 
parison to our models at [Fe/H] = -2.3. The SAGA database 
contains 287 stars with C measurements in this range, out of 
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which 59 are CEMP stars. This amounts to a CEMP fractior0 
of 20.6%. Adding the additional CEMP stars from Masseron' 
20T0| and |Lucatello et al.| ( |200 6 ) yields 92 CEMP stars 



et al. 



and 2 NEMP stars, both of which are also (marginally) CNEMP. 
Despite a nominal /nemp of 0.022 in this range, we cannot ex- 
clude real fractions up to 0.088 with more than 99% confidence. 

The NEMP/CEMP ratio appears to be larger at both lower 
and higher metallicity. The difference at lower metallicity is 
likely to be real, since at least 5 out of the 7 NEMP stars at 
[Fe/H] < -2.8 appear to be genuine. We can compute the 
probability that we observe at least 5 NEMP stars among 26 
CEMP-hNEMP stars given a real fraction / in a similar way, 
and find we can exclude / < 0.097 with 90% confidence and 
/ < 0.052 with 99% confidence. On the other hand, the status 
of the 3 NEMP stars at [Fe/H] > -1.8 is less certain and the 
high nominal NEMP/CEMP ratio could just be due to chance 
and small-number statistics. We can take these stars into account 
by combining them with the -2.8 < [Fe/H] < -1.8 sample, see 
the lower line of Table |2 We then find upper limits to /nemp that 
are only slightly larger than from the -2.8 < [Fe/H] < 1.8 sam- 
ple by itself, while the limits to /cnemp are even slightly lower. 

We conclude by noting that the SAGA database is by no 
means a statistically complete sample. Bright giants are clearly 
overrepresented with respect to fainter giants and turnoff stars, 
compared to what is expected from the relative lifetimes of 
these phases ( Izzard et al. 2009 ). However, neither the observed 
sample nor our model results show a strong dependence of the 
CEMP fraction on evolution state (represented e.g. by the effec- 
tive gravity). Therefore the results we describe below are proba- 
bly not greatly affected by this selection effect. 



3. Population synthesis models 

We have simulated populations of metal-poor halo stars in bi- 
nary systems using the rapid synthetic binary nucleosynthesis 
code of Izzard et al. ( 20041 "2006). The algorithm and the as- 



sumptions made in these simulations are discussed extensively 
in Paper I. We adopt a metallicity Z - 10 ([Fe/H] = -2.3) and 
select stars with ages between 10 and 13.7 Gyr (roughly corre- 
sponding to the age of the halo) and log g < 4.0. Our synthetic 
populations contain only binaries (see Sect.|5]for a discussion of 
the binary fraction). However, the predicted NEMP/CEMP frac- 
tion is independent of the binary fraction, assuming both CEMP 
and NEMP stars are formed through binary mass transfer. 

We consider models A, G and H from Paper I, to which we 
refer for details. Model A is our default model, where we param- 
eterize AGB evolution, including the efficiency of third dredge- 
up (TDU), according to the detailed models of Karakas et al.| 
( [20521 1. This yields C-rich AGB stars and thus CEMP progeni- 
tors for initial primary masses Mi > 1. 2 Mq. Hot bottom burning 
produces N-rich AGB stars for Mi > 2.7 M© in these models. 
In Model G we assume much more efficient TDU in low-mass 
AGB stars by reducing the minimum core mass and envelope 
mass required for dredge-up, as described in Paper I. This re- 
sults in almost all stars with initial masses 0.85-2.7 M© becom- 
ing C-rich AGB stars and thus potential CEMP progenitors. On 
the other hand, more massive stars - and thus NEMP progen- 
itors - are not affected by this choice. We note that more effi- 



' In 
a simi 



Izzard et al. 



2009 1 we found a larger CEMP fraction, 31%, from 



ar selection on [Fe/H] and log^. However, the dataset used in 
that paper combined the SAGA database with a selection of the data 
from [Lucatello et al.| ( |2006) which is strongly biased towards CEMP 
stars. 



cient TDU is also required to matc h the carbon-star lu minosity 
functions in the Magellanic Clouds (Izzar d et al.|2004 V Models 
A and G both assume efficient thermohaline mixing of the ac- 
creted matter throughout the stellar envelope of the secondary. 
Model H is similar to Model G but we do not allow thermoha- 
line mixing, such that accreted material remains at the surface of 
the secondary with unchanged abundances until first dredge-up 
occurs. 

In this study, in addition to varying the main assumptions 
about the input physics, we also vary the initial binary parame- 
ter distributions as detailed below. We have no direct informa- 
tion about these distributions in the halo at low metallicity, and 
must necessarily make some rather crude assumptions. Models 
1^ are alternative representations of the binary population in 
the solar neighbourhood, with which we test the sensitivity of 
our results to uncertainties in these distributions. 



1. 



Default model, as in Paper I; the initial primary masses Mi 
are distributed according to the solar neighbourhood IMF as 
derived by Kroupa et al. ( 1993[ l, the initial separations are 
drawn from a flat distribution in log a (with a between 3 and 
10^ Rq) and the initial mass ratios from a flat distribution in 
q - M2/M1. The IMF and mass-ratio distribution for this 
model are shown as solid lines in Fig. |2j5 and c. A flat dis- 
tribution in log a provides the best description of the binary 
population of young OB associations (Kouwenhoven et al. 
2007| l, and thus represents intermediate-mass binary popula- 
tions in the Galaxy (with Mi S 3 Mq). 
Primary masses and mass ratios are distributed as in Model 
1, but orbital periods are drawn from the distributions de- 
rived by|Duquennoy & Mayor (1991 hereafter DM91) for 



the local population of G dwarfs (with Mi ~ 1 Mq), i.e. 
a log-normal period distribution with a broad peak at 170 
years. 

As model 2, but here also the mass ratios are drawn from the 
distribution proposed by DM91, i.e. a Gaussian mass-ratio 
distribution with a broad peak at ^ = 0.23. This is shown as 
the dashed line in Fig.|2|;. 

Primary masses are taken from an alternative IMF for the 
solar neighbourhood by Miller & Scalo ( 1979 1, described by 
a log-normal distribution. 



dN 



d log Ml 



exp 



(log Ml -logMmed) 



(2) 



with a median mass of Mn,ed = 0.1 M© and a width ctm - 
0.67, while separations and mass ratios are distributed as in 
Model 1. 

Models 5 and 6 test the effect of two alternative forms of the IMF 
that have been suggested in the literature. 

5. As Model 4, but with a modified log-normal IMF (Eq. 2]i 
with Mtned - 0.79 M q and ctm - 0.5, as proposed by 
[LucateUo et al.| ( |2005a| l. This IMF is shown as a dashed line 
inFig.|2f). 

6. As Model 4 using a log-normal IMF with a much larger me- 
dian mass of 10 Mq and a narrower width, (Tm - 0.33, as 
proposed by Komiya et al. ( 2007| l. 



We also test the effect of varying the IMF parameters in 
Eq. (j2|i, Mtned and ctm, in a continuous manner between 0.1 Mq 
and 20 Mq and between 0.1 and 1.0, respectively. 
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Kom^O? : 




-0.5 0.5 1 1 3 



log Mj dN/dq 

Fig. 2. (a) Schematic regions in the parameter space of initial 
primary mass Mi and mass ratio q for the progenitors of EMP 
stars (M a; 0.8 M^, between the solid lines) and of CEMP and 
NEMP stars (light and dark shading, respectively); see text for 
details, (b) Three examples of the initial primary mass function 
( |Kroupa et al. 1993 , Lucatello et al. 2005a , Komiyaet al. 2007) . 
(c) Comparison of two initial mass-ratio distributions. 



4. Results 

Figure [2^ shows in a schematic way where the progenitors of 
CEMP and NEMP stars are located in the parameter space of 
initial primary mass M\ and initial mass ratio q. The constraints 
on age and gravity mean that we select stars with masses in a 
narrow interval around 0.8 M©. These stars are either primary 
components in binaries (indicated by the vertical lines in Fig.|2^) 
or the original secondaries in binaries where the more massive 
primary has evolved into a compact remnant. The latter have 
M2 - qM\ X 0.8 Mq, which defines a small range of mass ra- 
tios oc 1 /Ml for each Mi (shown by the curved lines^ The 
CEMP and NEMP progenitors lie in the grey-shaded regions in 
Fig. |2^, where the mass limits shown correspond to Model A. 
Only binaries with orbital periods in a certain range evolve into 
CEMP or NEMP stars; this is not shown, but our models indi- 
cate that this period range is not very sensitive to Mi . Another 
simplification in Fig.|2]is that we have ignored accretion: in re- 
ality the initial mass ratios of CEMP and NEMP progenitors are 
somewhat smaller because the secondaries accrete up to 0. 1 M© 
of primary material. 

This figure illustrates several effects that are visible in the re- 
sults we obtain. First, both the CEMP fraction among stars with 
M 0.8 Mq and the NEMP/CEMP ratio depend on the shape 
of the IMF, both increasing as the median mass of stars formed 



' Note that in Fig. [2] the indicated width AM of the mass interval 
around 0.8 Mq is arbitrary and only for illustration purposes; this can- 
cels out when calculating stellar number ratios as long as AM <k M. 



Table 3. Predicted number fractions of CEMP, NEMP and 
CNEMP stars among halo stars at [Fe/H] ~ -2.3 and log^ < 
4.0. The last two columns give the fraction /nemp of NEMP stars 
among stars that are either CEMP or NEMP and the fraction 
/cNEMP of CNEMP stars among CEMP stars. 



model 


CEMP 


NEMP 


CNEMP 


/nemp 


/cNEMP 


Al 


2.33 % 


0.27 % 


0.10% 


0.107 


0.041 


A2 


2.75 % 


0.34 % 


0.11% 


0.114 


0.041 


A3 


3.50% 


0.54 % 


0.18% 


0.139 


0.051 


A4 


3.18% 


0.48 % 


0.17% 


0.137 


0.051 


A5 


4.87 % 


1.04% 


0.36 % 


0.188 


0.074 


A6 


13.3 % 


19.7 % 


5.4% 


0.714 


0.406 


Gl 


9.5% 


0.27 % 


0.10% 


0.027 


0.011 


G2 


11.1 % 


0.34 % 


0.12% 


0.030 


0.011 


G3 


8.5% 


0.53 % 


0.19% 


0.060 


0.023 


G4 


10.6% 


0.48 % 


0.18% 


0.044 


0.017 


G5 


12.7 % 


1.04% 


0.39% 


0.078 


0.031 


G6 


15.6% 


19.7 % 


5.8% 


0.667 


0.372 


HI 


15.8% 


0.43 % 


0.30 % 


0.027 


0.019 


H2 


18.7% 


0.55 % 


0.38 % 


0.029 


0.021 


H3 


15.0% 


0.87 % 


0.60% 


0.057 


0.040 


H4 


17.9 % 


0.77 % 


0.53 % 


0.042 


0.030 


H5 


22.1 % 


1.68% 


1.16% 


0.074 


0.052 


H6 


43.5 % 


35.8 % 


22.8 % 


0.634 


0.525 



increases. Second, since Aq oc 1/Mi, the fraction of binary sys- 
tems at a certain Mi that evolve into a CEMP or NEMP star 
decreases with primary mass. For a flat mass-ratio distribution 
this fraction is, to first order, inversely proportional to Mi. This 
effect results in a smaller NEMP/CEMP ratio than one would 
expect on the basis of the IMF alone. Third, NEMP progenitor 
systems have smaller initial mass ratios (O.I ^ q ^ 0.3) than 
CEMP progenitors (0.3 <, q £ 0.7 for Model A). Therefore the 
NEMP/CEMP ratio is sensitive to the initial mass-ratio distri- 
bution. On the other hand, because the orbital periods of CEMP 
and NEMP progenitors are similar, the NEMP/CEMP ratio is not 
sensitive to the initial period distribution. 

In Table |3] we give the number fractions of CEMP, NEMP 
and CNEMP stars resulting from our models, for various as- 
sumptions regarding the physical ingredients (models A, G and 
H) and the initial distributions of binary parameters (1-6). We 
also compute the (C)NEMP/CEMP fractions /nemp and /cnemp 
from these models for comparison with the observational con- 
straints derived in Sect. |2.1| Since our models are for [Fe/H] = 
-2.3, this comparison should strictly be made only with the sub- 
sample having -2.8 < [Fe/H] < -1.8. A substantial fraction 
of NEMP stars in our models are CNEMP, because even though 
intermediate-mass AGB stars convert most dredged-up carbon 
into nitrogen by HBB, these stars still produce large amounts 
of carbon. This material can be further diluted by thermohaline 
mixing in the secondary and therefore in model sets A and G a 
smaller proportion of NEMP stars are CNEMP than in model set 
H, in which no such dilution occurs until first dredge-up. 

The results for models with default binary parameter distri- 
butions (Al, Gl and HI) are discussed extensively in Paper I. 
Our default model (Al) fails to account for the large observed 
CEMP frequency, while Model Gl yields an increase of the 
CEMP fraction to almost 10%, substantially closer to the ob- 
served value. This coiTesponds to decreasing the lower primary 
mass limit of CEMP progenitors in Fig. |2] to about 0.85 Mg. 
Inhibiting thermohaline mixing increases the CEMP frequency 
to more than 15% (Model HI). In all these default models the 
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model A 



model G 



CEMP fraction 



ii4 



)K5 



> 10% 

> 15% 

> 20% 

> 25% 




< 0.116 (99.9%) 
NEMP/CEMP < 0,088 (99.0%) 

< 0.057 (90.0%) 



1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



Fig. 3. Dependence of the CEMP fraction (upper panel) and the 
fraction /nemp of NEMP stars among all CEMP+NEMP stars 
(lower panel) for our default model A, as a function of the 
median mass M^ed and dispersion cr^ for a log-normal IMF. 
Contours in the upper panel correspond to CEMP fractions be- 
tween 10% and 25%, which spans the range of observed values. 
The soUd contours in the lower panel correspond to the 90%, 
99% and 99.9% confidence upper limits on /nemp in the range 
-2.8 < [Fe/H] < -1.8 and logg < 4, while the dashed line 
shows the 99% upper limit on /cnemp (see Ta ble [2]l. The star 




symbols labelled 4, 5 and 6 indicate the IMFs of 
1979| l, |Lucatello et al.| ( [2005a| l and |Komiyaet a 



Miller & Scalo 



(20071. 



small overall NEMP fraction (< 0.5%) is consistent with the ob- 
servations, and Models Gl and HI with enhanced dredge-up can 
comfortably reproduce the observed small /nemp- Based on the 
considerations in Sect. 2. 1 Model Al can be rejected with more 
than 99% confidence on the basis of its /nemp value - but this 
is entirely due to the small CEMP fraction in this model. If we 
consider only CNEMP stars, the constraints are weaker and also 
Model Al (with /cnemp = 0.041) cannot be excluded with more 
than 90% confidence. 

4. 1 . Variations in parameter distributions of tiie local binary 
population 

Distributions 2 and 3 give some insight into the dependence of 
our results on uncertainties in the local (solar-neighbourhood) 
binary parameter distributions. Model A2 shows a modest in- 
crease in both the number of CEMP stars and NEMP stars com- 
pared to Model A 1 , because the peak in the |Duquennoy & Mayor| 
( I99I|l period distribution coincides with the period range in 
which mass transfer is most effective. This is partly compen- 
sated by the wider range of periods in the DM91 distribution. 




Fig. 4. Same as Fig.^for model G, with enhanced dredge-up in 
low-mass AGB stars. 



Similarly modest effects are seen when considering Models G2 
and H2. /nemp and /cnemp are hardly affected by this choice. 

Adopting also the DM91 mass-ratio distribution in 
Model A3 gives a further increase in the CEMP and NEMP frac- 
tions, by factors of 1.5 and 2 compared to Model Al, as this dis- 
tribution favours the initial mass ratios with which most CEMP 
stars and all NEMP stars are born in model set A (see Fig.|2j;). A 
similar increase in the number of NEMP stars is seen in model 
sets G and H. On the other hand, in Model G3 the number of 
CEMP stars is somewhat smaller than in Model Gl because most 
CEMP stars now come from binaries with primary masses 0.85- 
1.2 Mq and q > 0.7, which is disfavoured by the DM91 mass- 
ratio distribution (Fig. |2};). The same applies to Model H3. In 
Models G3 and H3 /nemp and /cnemp increase by a factor of 
2 compared to the default distribution. While less likely to be 
correct, these models cannot be excluded with more than 90% 
confidence (Table|2]i. 

Adoption of the Miller & Scalo ( 1979) IMF (Models A4, G4 
and H4) gives somewhat higher CEMP and NEMP fractions than 
the |Kroupa et al.| ([T993) IMF, but the effect is modest: less than 
factors of 1.4 and 1.8 for the CEMP and NEMP stars, respec- 
tively. We conclude that uncertainties in the distributions of the 
solar-neighbourhood binary populations affect the CEMP and 
NEMP fraction and their ratio by at most a factor of 2. 



4.2. Modified forms of the initial mass function 

Our results for distributions 5 and 6 show the effect of a mod- 
ified initial mass function. Model A5 assumes the IMF sug- 



gested by Lucatello et al. (2005a) as being required to reproduce 
the large CEMP fraction. The fraction of CEMP stars increases 
by a factor of 2 compared to the KTG93 IMF (Model Al), as 
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model H 




Fig. 5. Same as Fig. |3]for model H, without thermohaline mix- 
ing. 



this IMF favours intermediate-mass stars over low-mass stars 
(Fig. However, the CEMP fraction still falls sh ort of the 
observed value. Our CEMP fraction is smaller than ILucatellol 
|et al.| ( |2005a| l predict because our ranges of initial primary mass 
and orbital period that produce CEMP stars are more restrictive. 
With enhanced third dredge-up the increase in the number of 
CEMP stars is only a modest factor of 1.3 (Model G5) to 1.4 
(Model H5). Models A5-H5 also show an increased (C)NEMP 
fraction, by a factor of 4 over models Al-Hl, respectively, 
which is the direct result of a larger weight of intermediate-mass 
stars in this IMF Models G5 and H5 have /nemp ~ 0.07-0.08 
and /cNEMP ~ 0.03-0.05, which is still (marginally) compatible 
with the observed limits and cannot be excluded with high confi- 
dence. We note that model H5, which combines efficient dredge- 
up, inhibited thermohaline mixing and a modest change to the 
IMF, is the only model that comfortably reproduces a CEMP 
fraction of about 20% while remaining consistent with the small 
observed number of NEMP stars at [Fe/H] > -2.8. 

The effects noted above are much more extreme when we 
assume the IMF suggested by [Komiya et al. ( 2007 [ l, which has 
a median mass of 10 Mq. This IMF gives rise to a substantial 
CEMP fraction, even with default input physics (Model A6), al- 
though it remains smaller than 20% even in Model G6: the effect 
of enhanced dredge-up is only modest in combination with this 
IMF. However, the main effect of such a top-heavy IMF is an 
enormous increase in the number of NEMP stars. One would 
expect about 20% of metal-poor stars to be NEMP and more 
than 5% to be CNEMP according to Models A6 and G6, and 
these numbers are even higher (36% NEMP and 22% CNEMP) 
in Model H6. In these models /nemp and /cnemp are strongly in- 
compatible with the observed limits, at least for [Fe/H] > -2.8. 



On the other hand, given the much larger apparent /nemp 
among stars with [Fe/H] < -2.8 one might more plausibly con- 
sider a top-heavy IMF at these low metallicities. Although we 
cannot really compare our models to this low-metallicity sam- 
ple, we might assume that the results in Table[3]are not sensitive 
to the metallicity and also apply to [Fe/H] < -2.8. If this were 
the case, then using the analysis of Sect. |2.1| we find the prob- 
ability of finding at most 7 NEMP stars among the sample of 
26 CEMP and NEMP stars to be only 3 x 10 ^ 4 x 10"^ and 
2 X 10 for Models A6, G6 and H6, respectively. Considering 
/cnemp, the probability of finding 5 CNEMP stars among the 24 
CEMP stars would be 3.5%, 7% and 0.16% for Models A6, G6 
and H6, respectively. This crude analysis suggests that the data 
at [Fe/H] < -2.8 do not support a Komiya et al. (2007i top- 
heavy IMF, unless many non-C-rich NEMP stars have remained 
undiscovered. 

We investigate the constraints on the IMF further in Figs. [3]- 
|5] where we convolve our model calculations for Models A, G 
and H, respectively, with a log-normal IMF of the form (j2|i using 
the median mass and width of the distribution as free parame- 
ters. The top panels show contours of the CEMP fraction while 
the bottom panels show contours of /nemp, both as a function 
of Mmed and ctm- The contours in the lower panels correspond 
to the 90%, 99% and 99.9% confidence limits on /nemp and the 
99% limit on /cnemp (dashed line) for -2.8 < [ Fe/H] < -1.8, 
see Table |2] This is analogous to the analysis of |Komiya et al.| 
( 2007[ l, the main difference being that we use our detailed pop- 
ulation synthesis models as input and that we also consider the 
NEMP/CEMP fraction /nemp as a consti-aint. 

Fig. [3] demonstrates that Model A can only reproduce the 
observed CEMP fraction for M^ed > 1 Mq and small values of 
(Tm, i e. only for very narrow, top-heavy forms of the IMF. On 
the other hand, the upper limit on /nemp only allows small Mmed 
and again favours small values of cr^. In order to reproduce both 
with our default model, we would require a highly fine-tuned 
IMF with Mmed ~ 2M0 and ctm ^ 0.2, i.e. one which only 
produces intermediate-mass stars in a small mass range. 

These stringent constraints are relieved with Model G, which 
produces many more CEMP stars but essentially the same num- 
ber of NEMP stars, and thus a lower /nemp- In this case a wider 
range of IMFs with Mmed < 2 M© is able to reproduce /nemp, 
although a large CEMP fraction (> 15%) still requires a rather 
top-heavy IMF. In this case, however, CEMP fractions > 10% 
are also produced by a variety of broad, bottom-heavy IMFs. 
There is a band of values of M^ed and cr^, which includes the 
Miller & Scalo and Lucatello IMFs, that allows CEMP frac- 
tions > 10% with acceptable /nemp values. In Model H both the 
CEMP and NEMP fractions increase with respect to Model G, 
while /nemp remains almost unchanged. This model allows a 
range of bottom-heavy IMFs (M^ed < 2 Mq) that reproduce both 
/nemp and a CEMP fraction of at least 15%. If we consider only 
CNEMP stars (dashed lines in Figs. [3]-^ the ranges of allowed 
IMF shapes do not change substantially, except that for model G 
a range of broad IMFs with M^ed > 2 M© cannot be firmly ex- 
cluded. 

We conclude that only bottom-heavy IMFs are potentially 
able to reproduce both the observed CEMP fraction and the 
NEMP/CEMP fraction among stai's with -2.8 < [Fe/H] < -1.8, 
provided that efficient dredge-up of carbon occurs in low-mass 
AGB stars. On the other hand, top-heavy IMFs (with Mmed > 
2 Mq) are unable to produce a low enough /nemp, if indeed HBB 
occurs in metal-poor AGB stars with M > 2.7 Mq as detailed 
models indicate. 
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5. Discussion 

In the following subsections we discuss how various uncertain- 
ties in our models and assumptions could affect our results. 



5.1. The lower mass limit for hot bottom burning 

The likelihood of NEMP star formation depends on which 
masses of star undergo hot bottom burning. The higher the min- 
imum initial mass for the onset of HBB, the fewer NEMP stars 
will form. Our results are based on an adopted lower mass limit 
for HBB, Mhbb, of 2.7 Mq in accordance with the results of 
[Karakas & Lattanz iol (|2007|. On the other hand, ISuda e t aT] 
(2011) argue that if Mhbb is larger than 5 Mq then the observed 
NEMP/CEMP ratio for [Fe/H] < -2.5 is compatible with the 



low metallicity it seems likely that some extra mixing process is 
at work. (|Stancliffe et al.|2009tpa"sseron et al.|2010[|Lucatello| 



top-heavy IMF of |Komiya et al.| ( |2007 1. 

The lower mass limit for HBB is thus an important quantity. 
Unfortunately, there is some uncertainty in this quantity, with 
differ ent stellar evol ution codes giving different values. At Z - 
10"'^,|Herwigl(|2004| finds Mhbb ~ 3.5 Mq for his models which 



include convective overshooting. A similar value is obtained by 
Stanclifi'e & Glebbeek](j20Q8 ), though they find some CN cycling 



occurring in a 3 Mq model. Karakas (201 0) finds mild HBB to 
occur in a 3 Mq model at this metallicit y. At a higher metallicity 



of Z = IO~^, |Weiss & Ferguson! ( |2009| l find hot bottom burning 
does not set in until the initial mass is > 5 Mq. 

Several groups have made grids of AGB models over a range 
of masses and metallicities. It is instructive to compare their re- 
sults, as this gives some indication of the uncertainties in the 
determination of Mhbb- Campbell & Lattanzio] ( |2008j l find that 
HBB occurs for masses equal to 2Mo and above for all the 
metallicities they study: from [Fe/H] = -3 right down to the 
metal-free population. Lau et al. (2009i find Mhbb ^ 3Mo for 
models without overshooting for metallicities from [Fe/H] - 
-2.3 down to [Fe/H] = -6.3, the lowest metallicity in their 
grid. When they include convective overshooting during the pre- 
AGB evolution Mhbb drops to 2Mo. In stark contrast to this 
is the work of Suda & Fujimoto ( 2010| l who do not find HBB 



in any of their grid of models, which spans metallicities from 
[Fe/H] = -2 down to [Fe/H] - -5 (plus Z - 0) and masses 
from 0.8 Mo to 9 Mq. 

The origin of these differences is not easy to determine. 
Certainly one contributing factor is that different groups use dif- 
ferent mixing lengths. This will change the efficiency of convec- 
tion and hence the mass at which HBB sets in. It seems likely 
that the efficiency of third dredge-up (one of the perennial uncer- 
tainties of AGB modelling) also affects the occurrence of HBB, 
with more efficient TDU likely to cause the onset of HBB at 
lower masses. Another possibility is the treatment of molecular 
opacities. Marigo (2007 ) found that when a star reaches C/O > 1 
then hot bottom burning can be quenched in some cases. This 
would shift Mhbb to higher masses, although |Marigo| ( [2007[ ) ex- 
pects the effect to be small for Z < 10"^. 

5.2. Non-canonical evolution 

The preceding discussion has focused solely on canonical mod- 
els of stellar evolution, that is to say, models that only trans- 
port material by convection. But it has long been thought that 
some extra mixing process may be at work in AGB stars, circu- 
lating material from the base of the convective envelope to the 
hydrogen burning shell. Whether this is necessary or not is an 



|et al.|201 1|). The physical cause of this extra mixing is unknown. 
Stancliffe (2010) has shown that thermohaline mixing (as cali- 
brated to match the abundance changes on the RGB) is unable to 
affect nitrogen production in low-metallicity AGB stars. 

The effect that extra mixing would have on nitrogen pro- 
duction is difficult to assess while the nature of the mechanism 
remains unknown. From observations of CEMP stars it seems 
that only a weak circulation is required, as these stars still show 
carbon enhancements in excess of nitrogen enhancements, and 
few CEMP stars show '^C/'^C ratios as low as the CNO cycle 
equilibrium value. So while extra mixing appears to be needed, 
it seems unlikely that it will significantly change the minimum 
mass for HBB. 

At low metallicity, AGB models predict mixing episodes not 
seen at higher metallicity. The intershell convection zone is able 
to penetrate into the H-rich envelope, sucking protons down to 
high temperatures and leading to a large release of energy from 
hydrogen burning. These dual shell flash events (or proton inges- 
tion episodes) can lead to the intershell convection zone splitting 
into two regions. One is driven by helium burning and the other 
by hydrogen burning. In the hydrogen-burning driven convec- 
tion zone, the CN cycle takes place and potentially the nitrogen 
abundance can be elevated when the convective envelope pen- 
etrates into this region. Several groups have recently attempted 



to model this phase, including Campbell & Lattanzio (2008|, 
Lau et al.' (20^, 'Iwamoto| ( .2009) , .Cristallo et al.| ( ,2009 ) and 
Suda & Fujimoto (20l0|l. However, these models are all based 
on ID mixing length theory which is almost certainly incorrect 
in these circumstances. What is really needed are hydrodynam- 



ical simulations of this phase of evolution. Recently, |StancUffe 
et al. ( 201 I| l performed such simulations. They found serious 



open question at higher metallicities (see [Karakas et al. 
and [Busso et al.||2010] for both sides of the argument). 



2010 



mt at 



deficiencies in the way ID evolution codes model the process of 
proton ingestion, noting that the transport of hydrogen is not dif- 
fusive in character, and that the velocities of the flows involved 
are significantly in excess of those given by mixing length the- 
ory. Crucially, they also found no evidence of the convective re- 
gion splitting into two zones, despite significant energy genera- 
tion from hydrogen burning. 

While dual shell flash events may enhance the nitrogen abun- 
dance, it seems unlikely that they would produce nitrogen-rich 
AGB stars, as this would require a large degree of hydrogen 
ingestion together with fairly complete CN cycling. Whether 
this is the case will have to remain unknown until more de- 
tailed multi-dimensional calculations can be made. In summary, 
it seems improbable that non-canonical evolution will produce 
nitrogen-rich AGB stars over a greater mass range than currently 
predicted. 

5.3. Mass transfer efficiency 

In our models mass transfer is treated as a combination of 
Roche-lobe overflow (RLOF) and Bondi-Hoyle wind accretion 
( [Bondi & Hoyle||1944| l. The modelling of both processes has 
substantial shortcomings when applied to binaries containing 
AGB mass donors. Roche-lobe overflow is relevant in binary 
systems that are close enough for the AGB star to fill its Roche 
lobe, and in most cases leads to unstable mass transfer and a 
common envelope. This evolution channel produces only a small 
fraction of the CEMP and NEMP stars in our models. However, 
the classical Roche geometry that we have assumed may not ap- 
ply to binaries containing luminous AGB stars, where radiation 
pressure may strongly modify the sizes and shapes of the Roche 
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lobes ( |Schuerman|1972[|Dermine et al.|2009l l. This could affect 
the stability of RLOF and the number of CEMP stars and NEMP 
stars produced by the RLOF channel. Because AGB stars under- 
going HBB are intrinsically brighter than lower-mass AGB stars, 
the effect may be largest on the formation of NEMP stars. 

The dominant formation mechanism of CEMP and NEMP 
stars in our models is wind accretion. By adopting Bondi-Hoyle 
accretion we make the implicit assumption that the wind velocity 
exceeds the orbital velocity. This assumption breaks down when 
the mass-losing star is an AGB star in a relatively close orbit. 
AGB mass loss is driven by a combination of stellar pulsation 
and radiation pressure on dust. The wind is accelerated in the re- 
gion where dust condenses, at a distance of several stellar radii. 
Hydrodynamical simulations show that when the companion or- 
bits at a distance where the AGB wind has not yet been acceler- 
ated to its terminal velocity, the gas flow geometry resembles 



Roche-lobe overflow rather than a Bondi-Hoyle flow (Nagae 


etal.|20041 


de Val-Borro et al.|2009l|Mohamed & Podsiadlowski 


20071 120 11 


1. In this transition regime of 'wind Roche-lobe over- 



find that the companion can accrete up to 50% of the mass lost by 
an AGB star, a much higher accretion efficiency than expected 
from the Bondi-Hoyle mechanism (:$ 10%). Abate et al. ( 2012| l 
have applied a simple model for wind-RLOF in a population 
synthesis study and find that this mechanism increases the oc- 
currence of CEMP stars by a factor of up to 1.5 compared to 
Bondi-Hoyle accretion. 

Of particular relevance for this study is the question whether 
a more realistic treatment of mass transfer would change the 
NEMP/CEMP ratio, i.e. whether it would affect the formation 
of NEMP stars in a different way than that of CEMP stars. This 
question is difficult to answer without dedicated mass transfer 
calculations for these cases. Such calculation are difficult to per- 
form because mass loss of low-metallicity AGB stars is poorly 
understood and is not constrained by observations. If radiation 
pressure on dust plays a similar role as at solar metallicity, then 
one might expect a difference between massive AGB stars that 
undergo HBB and lower-mass AGB stars that do not. Low-mass 
AGB stars produce carbon-rich dust, which condenses easily at 
relatively high temperature and thus close to the star On the 
other hand, stars that have undergone HBB eject gas that is 
nitrogen-rich and poor in C and O, which may form dust much 
less efficiently and at much larger distance to the star, although 
little is known about dust formation under such circumstancefl 
Massive AGB stars may therefore accelerate their winds at larger 
distance than their lower-mass counterparts. In the framework of 
the wind-RLOF model, one may thus expect a larger increase in 
the number of NEMP stars than in the CEMP number If this 
presumption is correct, the observed NEMP/CEMP ratio would 
constrain any modifications to the IMF even more strongly than 
we find in Sect.|4] 

5.4. The binary fraction 

The results we present are based on an assumed binary frac- 
tion of 100%, i.e. we assume that the progenitor population of 
CEMP and NEMP stars consists only of binaries. For metal- 
poor halo stars very little is known with certainty about the mul- 
tiplicity fraction, but the predominance of binaries among the 
CEMP-s stars suggests that binaries are as common in the halo 
as they are among Population 1 stars in the Galactic disk. It is 
thus reasonable to base our discussion on the observed multi- 



An interesting test case is provided by 77 Carinae ( Gull 2008 1. 



plicity in Galactic disk populations, which is much better de- 
termined. Among nearby solar-type stars, IDuquennoy & Mayor 
(1991 ) found a multiplicity fraction of about 57% for solar-type 
stars, which was revised to 46 + 2% in a more recent analysis of 
a much larger sample by Raghavan et al.| ( 2010] ). These studies 
concern stars of somewhat lower mass than the progenitors of 
CEMP stars. Perhaps more relevant, certainly for NEMP stars, 
is the binary fraction among intermediate-mass stars which was 
studied by Kouwenhoven et al. ( 2007 1 in the young OB asso- 
ciation Sco OB2. They found a binary fraction of at least 70% 
(3-cr lower limit) and probably much closer to 100%. A similarly 
high overall binary fraction is found among massive 0-type stars 
(Mason et al.„2009; Sana et al. 2012j|. A comparison of these 
results suggests that the binary fraction among Galactic disk 
stars is an increasing function of stellar mass, a conclusion also 
reached by Lada (2006) for low-mass stars. If a similar trend ex- 
ists for metal-poor halo stars, then the predicted NEMP/CEMP 
number ratio from our models would be higher. This would im- 
pose an even stronger constraint on the IMF. 

6. Conclusions 

Our results show that the occurrence of NEMP stars, in partic- 
ular the observed ratio of NEMP to CEMP stars, sets impor- 
tant constraints on the initial mass function of the early Galactic 
halo. In the mass-transfer scenario, the production of nitrogen 
by hot bottom burning in intermediate-mass AGB stars implies 
that any shift in the IMF towards more massive stars is accom- 
panied by an increase in the NEMP/CEMP ratio. From the cur- 
rently known census of 12 candidate NEMP stars with \ogg < 4 
we derive an upper limit to the NEMP/CEMP ratio of 0.15 (at 
99% confidence). In the metallicity range -2.8 < [Fe/H] < 
-1.8 only two candidate NEMP stars are known, and we find 
NEMP/CEMP < 0.09 with 99% confidence. Our detailed AGB 
models at [Fe/H] - -2.3, in which HBB occurs in AGB stars 
with M > 2.7 Mq, exclude IMFs with median masses of ^ 2 M© 
for stars formed in this metallicity range. On the other hand the 
observed NEMP/CEMP ratio is quite compatible with the solar- 
neighbourhood IMF, provided that low-metallicity AGB stars in 
the mass range 0.85-1.2 Mq are capable of dredging up carbon 
and turning a low-mass binary companion into a CEMP star, as 
we concluded in Paper I. 

The sample of known NEMP stars mostly have [Fe/H] < 
-2.8, which suggests that the NEMP/CEMP ratio is a strong 
function of metallicity. At [Fe/H] < -2.8 the observed NEMP/ 
CEMP ratio is approximately 0.27 (with a 99% confidence up- 
per limit of 0.51). From currently available model sets of AGB 
stars it is not clear how the mass limits for HBB behave at 
such low metallicities, and therefore it remains unclear exactly 
what the implications are for the IMF among the most metal- 
poor halo stars. Assuming HBB occurs for similar masses as 
at [Fe/H] = -2.3, the observed NEMP/CEMP ratio allows - 
and possibly requires - a shift in the IMF towards intermediate- 
mass stars. Nevertheless, a top-heavy IMF such as suggested by 
Komiya et al. (2007i remains firmly excluded, unless HBB es- 
sentially shuts off at the lowest metallicities ( Suda et al.|201 l| l. 
A larger census of NEMP stars would be very valuable to put 
firmer constraints on possible changes to the IMF at low metal- 
licity. 
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Appendix A: Comments on individual NEIVIP stars 



A.4. CS30314-067 



In this section we comment on the reported observations for our 
potential NEMP stars. If these stars are truly polluted with ma- 
terial from hot-bottom burning companions they are expected 
to be rich in nitrogen, sodium and magnesium (e.g. |Stancliffe] 
|2009) ). Further processing of C and N because of extra mixing 
may occur in the star if is evolved beyond first dredge up. 

Alpha elements such as calcium and titanium are enhanced 
by about 0.3-0.4 dex in normal Galactic halo stars. Some lithium 
enhancement is also expected because of hot-bottom burning 
although there are many other processes through which it may 



be made (Iwamoto 2009 StanclifFe 2010). Heavy s-process ele- 



ments such as barium and lead are not thought to be produced in 
intermediate-mass AGB stars so are not expected to be signifi- 
cantly enhanced in NEMP stars t| Lugaro et al.|2012 1 



The typical l-cr errors on abundance measurements [X/Fe] 
are 0.2-0.3 dex. 



A.1. CS22949-037 

The abundances of carbon and nitrogen in CS22949-037 are re- 
ported by .Norris eTall (|200 TT), IDep agne et~aL] p002] l, |Cayrel| 
[eTaL] ( |2004| ) and |Coheri etT l. (2008 ) who all find [N/Fe] > 1, 
[C/Fe] > 1 and [N/CI > 0.5. C S22949-037 is thus a CEMP and 
NEMP star. [Spite etal] ( |2006| l find ^^C/'^C = 4 which implies 
that the CN cycle is in equilibrium. These abundances are consis- 
tent with an intermediate mass AGB star undergoing hot bottom 
burning. The gravity of this star, log g - 1 .5 - 1 .7, implies that it 
has passed through first dredge-up and is at the point where ex- 
tra mixing begins but has not yet had time to significantly alter 
su rface abundances of c arbon and nitrogen (see Figs. 2, 6, and 7 
in lStanchff^e et al.|2009| l. 

Enhancements of sodium ([Na/Fe] - -1-1.57, 

-Hl.55, 



Andrievsky 



et al.| 2007 1 and magnesium ([Mg/Fe] = -1-1.55, lAndrievsky 
et al.||20l"0 i also suggest hot bottom burning is active. Barium 
is deficient ([Ba/Fe] = - 0.84, [Norris et al.|2001| -0.66, Cohen[ 



et al. 2008 -0.5, Andrievsky et al. 201 li and strontium is 



mildly enhanced ([Sr/ Fe] ^ -hO. 18, [Cohen et al.[|2008l h-0.17. 



[Andrievsky et al.[201 ip . These abundances fit the canonical view 
that most s-process enhancement occurs in low-mass AGB stars. 

A.2. CS22960-053 

This star is enriched in carbon and nitrogen with [C/Fe] - +2.05 
and [N/Fe] - +3.06 accor ding to |Aoki et"aL] ( |2007| l but in con- 
trast [John^reFal]([2007|) find [C/Fe] = [N/Fe] = h- 1.15. The 
discrepancy in nitrogen abundances may well be due to the re- 
spective use of CN and NH as abundance indicators. We can 
only conclude that this system may be a NEMP star The star 
is also enhanced in magnesium ([Mg/Fe] = 0.65) and barium 
([Ba/Fe] = 0.86). 

A.3. CS29528-041 

If we assume that the molecular abundance measurements of 
ISivarani et all ( [2006| , [CH/Fe] = 1.59, [CN/Fe] = 3.07 and 
[NH/Fe] = 3.00, are indicative of elemental abundances then 
this star satisfies our (C)NEMP criteria. With log g - 4.0 it is un- 
evolved and there is no mixing on the RGB. Sodium is enhanced, 
[Na/Fe] - 1.20, but magnesium is not, [Mg/Fe] = 0.40, espe- 
cially when compared to other alpha elements such as calcium 
and titanium which show similar enhancements. CS29528-041 
is rich in lithium, log eu = 1.71, and barium [Ba/Fe] = 0.97. 



CS303 14-067 classifies as a NEMP star according to the abun- 

[C/Fe] = -hO.5 and 



(2002 



(200711 instead find [C/Fe] = 



dances determined by Aoki et al. 
[N/Fe] = -Hl.20, but [Johnson"etaL] i 
-1-0.25 and [N/Fe] = H-0.50. Some of the discrepancy may result 
from the different nitrogen abundance indicators used and the 
different adopted gravities, although the [Aoki et al.[p002) study 
also has greater resolution. 



A.5. CS30322-023 



Masseron et al. (2006 1 have suggested that CS30322-023 is an 



AGB star because its gravity is low, \ogg = -0.3, but this is 
contradicted by [Aoki et aE] (|2007[) who claim log g = 1 .0. In 



either case extra mixing on the RGB may have occurred. Both 
[Masseron efa T (2006) and ' Aoki et"aL] ([2007^ claim [C/Fe] = 
-1-0.6 and significant nitrogen enrichment ([N/Fe] = -1-2.81 and 
-1-2.47 respectively) which suggest this is a NEMP star. Its abun- 
dances of sodium, magnesium and heavy elements are also en- 
hanced ([Ba/Fe] ^ -hO.6 and [Pb/Fe] = -i-1.49). 



A.6. HD25329 



HD25329 is an unevolved star {\ogg 
[N/Fe] = H-1.0 and '^c/'^c > 40 ([Fulbright 



4.6) with rC/Fe] = 0.1 



2000 



Gratton 



l et al.|2000[ and [Gratton et al.|2003| l. Its only sUght enhancements 
in sodium and magnesium ([Na/Fe] - H-0.24 and [Mg/Fe] = 
-1-0.59; Gratton et al. 2 003 ) suggest it has not been polluted by 
a hot bottom burning AGB star. The origin of nitrogen in this 
object remains unexplained. 



A.7. HD206983 

The metalUcity of HD 206983 is rather large, [Fe/H] = -0.99 
(Masseron et al.'2010), suggesting it is a nitrogen-rich CH-star 
rather than a NEMP star. This is supported by an enhancement 
of barium, [Ba/Fe] = h-0.92 ( [Masseron et al.|2010| l. 



A.8. HE0400-2030 

The carbon- and nitrogen-rich star HE0400-2030, with [C/Fe] = 
-1-1.14 and [N/Fe] = H-2.75 according to high-resolution ob- 
servations {R = 50, 000) of Aoki et al. (2007]), qualifies as a 
CNEMP star. Lower-resolution observations (R = 20, 000) of 
[Lucatello et al.[ ( |2006l ) find less enhancement, [N/Fe] = +1.0 
and [C/Fe] = H-0.8. Sodium ([Na/Fe] = -h0.71), magnesium 
([Mg/Fe] = -hO.62) and barium ([Ba/Fe] = -hI.64) are also en- 
hanced. 

If this star accreted from a hot bottom burning companion, 
then the accreted material must have been diluted during the 
main sequence (e.g. by thermohaline mixing, Stanchffe et al. 
[2007 ), as the sodium and magnesium abundances are lower than 
predicted by theoretical models. 



A.9. HE1031-0020 



[Cohen et al.[ ( |2006| find [C/Fe] = -h1.63 and [N/Fe] = -h2.48 
for HE103 1-002, i.e. it is a CNEMP star With logg = 2.2 
it is a moderately evolved giant and should have finished first 
dredge up without significant further extra mixing. Magnesium 
is slightly enhanced [Mg/Fe] - +0.5 but sodium is not mea- 
sured. Calcium, heavy s-process elements and lead are enhanced. 
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[Ca/Fe] = +1.12, [hs/ls] > -i-l and [Pb/Fe] = +2.66, suggest an 
odd evolutionary history for this star which does not necessarily 
involve a hot bottom burning companion. 

A 70. HE1 337+001 2 

The most extensive set of abundance measurements for the 
NEMP st ar HE1337-H0012 (also known as G64-12) is pro- 
vided by |Aoki et al.| p006) l who find [C/Fe] = -1-0.49 and 



[N/Fe] = H-1.42. The sodium abundance is low, [Na/Fe] = -1.1, 
so HE1337-I-0012 probably did not accrete material from a hot 
bottom burning companion despite its NEMP status. 

A.11. HE1410+0213 

This star is in the samples of both Cohen et al. ( 2006)) and 



Masseron et al. (2010 1 but the authors disagree on its surface 
gravity with logg = 3.5 and 2.0 respectively. This gives incom- 
patible abundances of [C/Fe] = +2.33 and [N/Fe] = -h2.94 



according to [Masseron et al. (20101 and [C/Fe] = -1-1.83 and 



[N/Fe] = +1.73 from [Cohen et al] ( |2006j ). Given the similar 
resolution of the studies it is unclear which is correct. 

A.12. HE1413-1954 



Lucatello et al. ( 2006[ l find this is an unevolved (logg > 4) 



CNEMP star with [N/Fe] = -h2.5 and [C/Fe] = -Hl.7. [Barklern| 
[eTaL] ( |2005] l measure logg = 4.59, [C/Fe] = -hI.45 and 
[Sr/Fe] = -0.47. The lack of further abundance measurements 
prevents a reliable determination of the origin of the nitrogen in 
this star. 



A. 13. HE21 50-0825 and HE2253-421 7 

HE2 150-0825 and HE2253-4217 are NEMP stars according to 



the carbon and nitrogen abundances measured by Lucatello et al. 
( |2006j who find [C/Fe] = +0.3 and [N/Fe] = +1.3, and 
[C/Fe] = +1.0 and [N/Fe] = +1.5 respectively. The status of 
HE2253-4217 as a NEMP is marginal with [N/C] = 0.5, while 
HE2150-0825 has a low total CN abundance, [(C +N)/Fe] = 
0.76. 

A. 14. CS29528-028 and SDSS1 707+58 

Neither CS29528-028 nor SDSS 1707+58 has measured nitro- 
gen, so they are not classed as NEMP stars, but their sodium 
and magnesium abundances are unusual. CS29528-028 is an un- 
evolved star with log g = 4.0. Aoki et aT] ([2007 ) report that 



CS29528-028 has [Na/Fe] = +2.68 and [Mg/Fe] = +1.69 de- 
spite a normal calcium abundance ([Ca/Fe] = +0.46). Its barium 
abundance is [Ba/Fe] = +3.27. 

SDSS 1707+58 is similarly enhanced with [Na/Fe] = +2.71 



and [Mg/Fe ] = +1.13 as well as [Ba/Fe] = +3.40 ( [AoE 
[et al.[[2008]l. This sta r was recently found to be an RR Lyr 



star by Kinman et al. ( 2012, ), who also find a lower metallic- 



ity [Fe/H] = -2.92 and somewhat lower sodium and barium 
abundances. 

Based on their sodium and magnesium abundances, these 
stars are quite possibly polluted by a companion undergoing hot 
bottom burning, although their barium enhancements do not sup- 
port this view. Unfortunately, until their nitrogen abundance and 
carbon isotopic ratios are measured the status of these stars re- 
mains unknown. 



